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Abstract. We present the first experimental confirmation of the so-called "self- 
phase-locked delay interferometry" . This laser frequency stabilization technique 
consists basically in comparing the prompt laser signal with a delayed version of itself 
that has been reflected in another LISA satellite 5 • 10^ m away. In our table-top 
experiment, the phase of a voltage controlled oscillator is stabilized by means of a 
control loop based on this technique. In agreement with the theory, the measured 
unity gain frequency is not limited by the inverse of the used delay (1.6 /is). In the 
time domain the system also behaves as predicted, including the appearance of a 
quasi-periodic "ringing" just after the lock acquisition, which decays exponentially. Its 
initial amplitude is smaller when the loop gain is slowly ramped up instead of suddenly 
switched on. 

E-mail: antonio . garciaSaei . mpg . de 

1. Motivation 

LISA is an ESA-NASA project to detect gravitational waves, involving three spacecraft 
flying in an equilateral triangle formation approximately 5 million kilometres apart. 
Together, they will act as a Michelson interferometer, covering a frequency range from 
0.1 mHz to IHz and having a typical strain sensitivity of 10~^^ pp. 

The phase noise of the LISA lasers would limit the sensitivity of the interferometer 
despite the use of traditional frequency stabilization techniques such as Pound-Drever- 
Hall with a stable reference cavity. TDI t2| represents an option to overcome this problem 
by postprocessing the acquired data. Nevertheless, it requires high performance of the 
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pre-stabilization methods mentioned before and it becomes more complicated when 
spacecraft motions are taken into account. 

A traditional approach used with ground based detectors consists in locking the 
laser frequency to the arms of the interferometer. The LISA arms are good candidates 
for this technique due to their exceptional stability in the measurement frequency band, 
but the delay caused by the roundtrip travel time between two satellites (33 s) had long 
been considered an insurmountable limitation. The control bandwith of this kind of 
loop is typically reduced to frequencies well below the inverse of this delay |3j , but that 
would mean a fraction of 30 mHz control bandwidth for LISA whereas very high gain at 
these frequencies is necessary to make the stabilization useful. Recently, some groups 
(see 13], [S] and jH]) have come up with control proposals and simulations achieving 
the necessary bandwidth and gain, what has been called the "self-phase-locked delay 
interferometry" . This paper describes an experimental demonstration of the principle 
of operation and the performance of the technique using an electrical model system. 

A voltage controlled oscillator (VCO) is stabilized in its frequency using a delay 
of r = 1.6 yus realized by 300 m of coaxial cable, and exhibiting the highest unity gain 
frequency (UGF) of the control loop beyond 1/r. 

The predicted noise suppression was confirmed by a direct measurement of the 
oscillator's signal. In the time domain, a quasi-periodic, exponential decaying transient 
that was predicted to appear just after the lock acquisition (JI] and f7|) could also be 
experimentally confirmed. Furthermore, its initial amplitude is reduced when the loop 
is closed by ramping up the gain instead of abruptly switching the loop on. 

2. Description and characterization of the system 
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Figure 1. Principle of operation of the self-phase- locked delay interferometry. 



Referring to Figure ^ an oscillator signal is split into two paths. After one of them 
has undergone a delay r, they are recombined, and a phasemeter detects their phase 
difference A$. When the loop is closed, A$ is used as the error signal r{t) for a control 
system (servo) with frequency response G{s). The output of the servo (feedback signal) 
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Table 1. Correspondence between the LISA properties relevant for the experiment 
and our prototype 





LISA 


Prototype 


Signal 


Laser 


VCO 


Delay r 


33 s 


1.6 /is 


1/r 


30niHz 


625 kHz 


freq. range 


0.1 mHz IHz 


2 kHz 20 MHz 



v{t) compensates the frequency fluctuations p(t) of the osciUator . Finally, q{t) is the 
frequency noise remaining in the stabilized system. In this paper the time will be called 
t and we will use the Laplace variable s = ic^;. Functions in the frequency domain are 
written with upper case letters and those in the time-domain in lower case letters. We 
will refer to the "system" as the whole stabilization loop, consisting of the "servo" and 
the "plant". 

In our experiment (see Figure El and Table P) the role of the LISA laser is played by 
a VCO working at approximately 72 MHz. Instead of the two times 5-10^ m pathlength 
between two LISA spacecraft, one of the signals goes through 300 m low-loss coaxial cable 
which causes a delay of 1.6 fis. The inverse of the delay is 625 kHz, and the frequency 
range equivalent to the LISA measurement window goes from 2 kHz to 20 MHz. 

In this section we will compare the theoretical response of the phasemeter A$ to 
frequency noise p{t) of the oscillator with the one measured in our prototype in the open 
loop case. Note that we actuate on the frequency of the oscillator instead of its phase, 
which results in an extra factor of 1// in the transfer function with respect to After 
that, we will discuss the open loop gain (OLG) and present the characteristics of the 
servo. 

2. 1 . Transfer function 

The transfer function of the plant without servo (see Figure CJ, measured from the 
frequency fluctuations of the oscillator p{t) (expressed in rad/s) to the phasemeter 
output A$ (expressed in rad) can be written as 

TT (■ N l-exp(-iwr) sin(a;r/2) . 

iftheo(it^) = : = T — -— exp(-iu;r/2). (1) 

Figure El shows the Nyquist and Bode representations of this function together 
with the data measured on the prototype. A model has been fltted to the data that 
additionally includes an extra delay r* of 75 ns. This delay r* accounts for effects in the 
VCO, the short interferometer arm and the phasemeter. The total effect is modelled by 
one single delay r* at the phasemeter output (see Figure E}- The model also includes 
additional poles at ^1,^3 and a zero at uj2 for the not ideal frequency response of the 
different components and is given by: 
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Figure 2. Transfer function of the plant. Left: Bode representation of tlie tlieoretical 
transfer function. The magnitude is given in units of r. Right: Nyquist representation. 
The curve labeled "theory" represents the theoretical transfer function. The curve 
labeled "data" represents the one measured on the prototype and the curve labeled 
"fit" shows the model presented in Equation El 



,. . sin(u;r/2) / • * 

HfiAiuj) =T— r^expl— iLur/ijexpl— icur , , . , , ^ , , , . , 

^ (W2) ' ' V + fJ\ ^^JV + fJ (2) 

with uji = 2tt- 530 kHz u;2 = 27r ■ 830 kHz cug = 27r ■ 12 MHz r* = 75 ns. 
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Figure 3. Table-top prototype: the phase measurement causes an extra delay r*. 




Figure El shows a more accurate description of the experimental setup including the 
extra delay r*. The transfer function shown in Figure |2l was measured as 

^ ^ y(iu;) + P(ia;) ^ ' 
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2.2. Open Loop Gain (OLG) 

The Open Loop Gain of the system was measured in a stable, well-behaved loop as 

The highest unity gain frequency (UGF) takes place at about 3.5 MHz (see Figure E}, 
clearly above the inverse of the delay (625 kHz). 




The noise suppression function predicted by such an OLG is given by: 

1 1 

^^"^^ " \l + OLG{iuj)\ " \l + G{iuo)H{iuj)\' 




For a better interpretation we consider the Nyquist representation of the OLG 
shown in Figure El The unity circle centred at (—1,0) represents the transition between 
noise suppression and noise enhancement. Besides, if the OLG encircles the (—1,0) 
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point, the system becomes unstable. In this representation, the additional phase lag 
r* causes a clockwise rotation of the spiral-like gain curve and thus limits both the 
loop bandwidth and the gain of the servo. It is important to distinguish between this 
spurious delay r* that appears in the control loop after the signals have recombined and 
the intended delay r that is applied to only one of the split signals. 

2.3. Controller 

Our controller (FigurelHl) consists of alternating poles and zeros at frequencies 106.6 kHz, 
172.4kHz, 843.7kHz, 1.8 MHz, producing a frequency response approaching /°'^ between 
200 kHz and IMHz. Besides, there is an extra integrator from DC to 100 kHz in order 
to increase the noise suppression in this region. 




Frequency [Hz] Frequency [Hz] 



Figure 6. Frequency response of the used servo and a fit to the measured data. 

The starting point for such a controller design is the theory presented in j3]. The 
only difference is that Sheard et al. consider the frequency actuator as part of the 
controller rather than of the plant. This corresponds to taking a factor of 1// from the 
transfer function presented here and giving it to the servo frequency response, which 
then becomes /~° ''. 

3. Results 

This section describes measurements of the noise behavior of the oscillator and compares 
them with the theory. The signal of the oscillator (72 MHz) is directly sampled at 
1 GHz. A software phasemeter performs a Single Bin Discrete Fourier Transformation 
on it (0,0) and delivers a time series of the oscillator's phase q{t) at 72 MHz data rate 
(see Figure El). We will study this signal in both the time and the frequency domain. 

3.1. Frequency domain 

For the measurement shown in Figure [71 white frequency noise {p{t) in Figure has 
been added into the system thus generating phase noise with a 1// linear spectral density 
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Figure 7. Linear spectral density (LSD) of the oscillator's phase (left) and frequency 
(right). The solid curves show the oscillator in "free-running" mode and the dotted 
ones refer to the stabilized state. The dashed curve on the right shows the frequency 
noise suppression predicted by the measured OLG ('Section |2.2|l . 



(LSD). Once the system is stabilized, the remaining noise Q{s) can be seen in the dotted 
curve which shows noise reduction at certain frequencies above 1/r. Superimposed to 
this curve, the disturbance sensitivity function 



is plotted, where P{iuj) is the LSD of the introduced frequency noise and OLG{iuj) is 
the measured OLG presented in Section ^ITA The reasonably good agreement between 
the two curves confirms the predicted noise suppression. 

3.2. Time domain investigations 

The time evolution of the oscillator's phase during lock acquisition is analyzed here. For 
all the figures of this subsection, the controller is turned on at t = and time units 
are scaled to r. For the solid curves, the gain of the controller was turned on abruptly, 
whereas for the dashed ones the gain was ramped up linearly during approximately 
16/is(10r). 

Figure IHl shows the transient for these two cases with white frequency noise added 
as described in Section EUl As predicted in |lj and [Zj, a pseudo-periodic transient can 
be observed just after the lock acquisition, whose initial amplitude is smaller when the 
gain is ramped up as opposed to the case of abrupt switching. 

The remaining noise is shown in the graph on the left in Figure IHl It does not show 
pure repetition but a structure typical to the sum of pseudo-harmonic narrow-band 
noise. This structure originates from the filtering done by the controller, as can be seen 
in Figure [7| On the right, 1// noise is injected instead of white noise. It shows how the 
system locks despite the strong perturbations that drive the phase of the oscillator over 
several hundreds of radians before the stabilization is turned on. 
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Figure 8. Lock acquisition as white frequency noise is being added into the system. 
The initial amphtude of the transient is smaUer when the gain is ramped up than when 
turned on abruptly. Left: General overview before and after the lock. Right: detailed 
view just after the lock. 
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Figure 9. Left: Detailed view of the time series shown in Figure |H1 It begins 200 r 
after locking. Right: Lock acquisition in the presence of 1// noise 



4. Discussion 

The "self-phase-locked delay interferometry" (see detects frequency fluctuations 
of a LISA laser by measuring the phase difference between the prompt laser signal 
and a delayed (r = 33 s) version of it that has been reflected on a different LISA 
satellite. In the hardware model presented here, the phase substraction takes place 
between the signal from a VCO and a second version of it, delyed by r = 1.6 /is. 
Frequency fluctuations of the VCO show up in our phase difference in the same way as 
frequency fluctuations of the laser do in the LISA configuration, which has allowed us 
the implementation of a frequency stabilization for the VCO based on the one described 
in Although it takes place in a different frequency range due to the small delay of 
1.6 /xs, it permits the experimental confirmation of the main features of the "self-phase- 
locked delay interferometry" . 

First of all, the highest UGF of this kind of stabilizations is not limited to values far 
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under 1/r, as was traditionally assumed (see [3 ). This can be seen from the measured 
OLG (Figure m and El) and from the noise plots (Figure [7|), in which frequency noise 
reduction takes place at frequencies above 625 kHz (1/t) and the highest UGF appears 
at about 3.5 MHz. 

Actually, the performance of our loop is only limited by the spurious delay 
T* = 75 ns present at the output of the phasemeter. This delay r* limits the bandwith 
and gain of the servo as discussed in section 12.21 Such delays appear frequently in the 
experimental realization of a phasemeter and therefore care should be taken to minimize 
them in further implementations of this technique. 

The frequency noise of the oscillator gets reduced when the stabilization is turned 
on, as can be seen in Figure [7| for the frequency domain and Figure |H1 and IHl for the time 
domain. This noise reduction is also in agreement with the performance that derives 
from the measured OLG (Figure Hj). 

Our hardware model of the "self-phase- locked delay interferometry" demonstrates 
that the 33 s delay present in LISA does not represent a fundamental limitation in the 
performance of the stabilization. The combination of the technique discussed here, 
together with TDI and its recent improvements ^01, may take us to a shot noise limited 
LISA without major hardware modifications in the actual baseline. 
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